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Abstract — Orthogonal Frequency Division Multiplex (OFDM) systems are 
very sensitive to frequency offset caused by tuning oscillator instabilities and 
Doppler shifts induced by the channel. The Schmidl and Cox Algorithm 
(SCA) is quite robust in estimating the frequency offset for systems with large 
OFDM symbol lengths. It uses two OFDM symbols for training with the first 
one having two identical halves. The frequency offset is estimated by correlat- 
ing a received sequence of samples equal to half the OFDM symbol length with 
the following received samples. The effect of Additive White Gaussian Noise 
(AWGN) in the estimation process is mitigated only if the number of samples 
used in the correlation, and hence the OFDM symbol size is large. However to 
be successfully applied to Broadband Fixed Wireless Access (BFWA) systems, 
OFDM should perform well even with smaller symbol lengths. In this paper 
we present the Residual Frequency Offset Correction Algorithm (RFOCA), 
which uses the SCA for initial frequency offset acquisition and follows it with 
a stage which reduces the initial residual frequency offset by tracking the 
phase of the decoded data. We show through simulation that the RFOCA 
achieves an error variance that is many orders of magnitude lower than at 
the end of the acquisition stage using the SCA alone. 



L Introduction 

Various solutions have been suggested as contenders to over- 
come tlie challenge of transmitting data to the subscriber at high 
data rates demanded by future broadband applications, for ex- 
ample Broadband Fixed Wireless Access (BFWA), Digital Sub- 
scriber Lines (DSL), Cable Modems and Satellite Communica- 
tion. This paper will address BfWA, owing to its advantages 
of ease of deployment and cost. Orthogonal Frequency Divi- 
sion Multiplex (OFDM) has been the most popular physical layer 
standard for BFWA systems due to its robustness and efficiency 
mainly because it uses a simple Fast Fourier Transform (FFT) 
of finite length A'' for modulation and demodulation. Among 
the standards that have placed their confidence in OFDM are 
HIPERLAN-2 [1] and IEEE 802.11-a [2]. 

The orthogonality of the consecutive OFDM symbols is main- 
tained by appending a length v cyclic prefix (CP) at the start of 
each symbol [3]. The CP is obtained by taking the last v samples 
of each symbol and consequently the total length of the transmit- 
ted OFDM symbols is {N -h v) samples. For each OFDM symbol 
to be independent and to avoid any Inter Symbol Interference (ISI) 
or Inter Carrier Interference (ICI), the length of the Channel Im- 
pulse response (CIR) should be less than xj -f 1 samples. Hence the 
distortion caused by the CIR only affects the samples within the 
CP. The receiver discards the CP and takes only the last N samples 
of each OFDM symbol for demodulation by the receiver FFT. In 
such a case, the linear convolution of the transmitted signal with 
the CIR is converted into a circular one. Consequently, for coher- 
ently modulated OFDM systems the effects of the CIR can then 
be equalized by an array of one- tap Frequency Domain Equalizers 
(FEQ) following the FFT. This is because the frequency selective 
fading channel can be approximated by a sum of flat fading chan- 
nels, provided the number of subchannels is large enough. Some- 
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Fig. 1 . Block Diagram of the OFDM System 



times a Time Domain Equalizer (TEQ) may be used to shorten 
the effect of the CER [4]. Figure 1 shows the block diagram of 
the system, where P/S and S/P represent the Parallel to Serial and 
Serial to Parallel functions respectively. 

Unfortunately, OFDM has been proven to be very sensitive to 
frequency offset and phase noise caused by tuning oscillator in- 
accuracies or Doppler shifts induced by the channel. The usual 
deployment scenario for a BfWA system employs a point to mul- 
tipoint configuration. Here, a single base station (BS) transmits 
data in short bursts to many subscriber units (SUs) placed at the 
user locations. The economic feasibility of a BFWA deployment 
depends heavily on the cost of the SUs. The use of low cost com- 
ponents, particularly the oscillators, is a major issue since their 
accuracy and stability are direcdy related to cost. Hence it is im- 
perative that if OFDM is to be used for BFWA transmission, it 
should be able to operate effectively using oscillators with only 
moderate performance and cost. 

An OFDM system effectively consists of N sinusoidal subcar- 
riers with frequency spacing 1/T, where T is the active sym- 
bol period of each subcarrier. The kth subcarrier will thus be at 
fk = /o + k/T, where /o is a reference frequency. Without loss 
of generality, we can assume /o — 0. The modulated subcarri- 
ers overlap spectrally, but since they are orthogonal over a symbol 
duration, they can be easily recovered as long as the channel does 
not destroy the orthogonality. An un windowed OFDM system has 
rectangular symbol shapes. Hence, in the frequency domain the 
individual subchannels will have the form of sine functions where 
the first sidelobe is only some 13 dB below the main lobe of the 
subcarrier. In the event of a frequency offset between the trans- 
mitter and the receiver the subchannels will not align. The effect 
of frequency offset is two fold, it will reduce the signal power of 
the desired subchannel outputs and also introduce ICI further in- 
creasing the interference levels. The ICI power can be significant 
even for small frequency offsets in OFDM due to the high sidelobe 
power of the subchannels. Thus the reduction of frequency offset 
is critical in OFDM systems. This paper is organised as follows, 



section II analyses the effect of frequency offset in OFDM analyt- 
ically. Section in briefly introduces the Schmidl and Cox Algo- 
rithm (SCA) [5] and section IV introduces the novel Residual Fre- 
quency Offset Correction Algorithm (RFOCA) that compliments 
the SCA and improves the performance. Section V shows results 
obtained by computer simulation and section VI concludes and 
proposes future work. 

II. Effect of Frequency Offset in OFDM 

All analysis and simulations in this paper are performed in the 
digital complex baseband domain. The nth sample of the mth 
OFDM symbol generated by the Inverse FFT (IFFT) at the trans- 
mitter is 



= /J E 0 < n < TV - 1 

fc=0 



(1) 



Am,k is the data symbol modulated on to the kih sub- 
carrier of the mth OFDM symbol. The data is con- 
verted into a serial sequence, then the CP of length v is 
added. Thus the mth transmitted OFDM symbol is s(m) = 
[sm,7v-i;,..,Sm,N-i,5m,o,",Sm,Ar-ir. The receiver discards 
the samples of the CP and uses the samples of the active OFDM 
symbol for decoding. We assume a finite length CIR with Nfi 
samples, h — [ho, /i7v^_i]^, where t; > AT^ - 1. We assume 
that A'' is large enough such that the frequency selective channel 
is divided into contiguous flat fading subchannels. If we assume 
that correct frame and timing synchronisation is achieved, then 
the received sequence after stripping the CP can be expressed as. 



^ fc=0 



(2) 



for 0 < n < A^ — 1. Here Hk is the transfer function of the 
channel at the subchannel index k,Wn is the Additive White Gaus- 
sian Noise (AWGN) term and e is the frequency offset relative to 
the intercarrier spacing. The symbol after FFT demodulation is 



N-l 



N 51 



(3) 



n=0 



Substituting (2) into (3) allows the output of the lib subchannel 
to be expressed as the sum of three components. 



sin(7rc) 
Nsm{7re/N) 



e^""'-^ +Ii-\-Wi (4) 



The first term on rhs of equation 4, denoted Ri hereafter, shows 
that the desired term Am.i^i experiences amplitude reduction and 
phase distortion. As A^ > ttc, N sinine/N) can be approximated 
by 7re. Hence the degradation of the required term vanishes when 
the received data does not experience any frequency offset. In 
practice, the effect of the channel transfer function Hi is removed 
by the FEQ following the FFT. Note that the FEQ is trained by 
transmitting a known training symbol occupying the OFDM sym- 
bol number Nt, A = [>ljVt,o, ^M.iv-i]^. Usually the training 
symbol is placed at the beginning of the frame. The ratio of the 
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Fig. 2, Analysis of CTIR vs Relative Frequency Offset, e 

decoded samples of the training symbol yivt,/* to the locally gen- 
erated copy of Ajsf^^i, where 0 < I < N — l/\s used to calculate 
the coefficients for the FEQ. 

The second term in equation (4) is the ICI and the third term is 
due to the AWGN. The ICI term is given as, 



- 2.(^-.'=^*Vsin7r(ifc + e-/)/Ar' 



fc=0 

fc/I 



(5) 



If we assume that the transmitted data is zero mean and uncor- 
rected (i.e. E[Am,k] = 0 and E[A^^iA*^J = \A\^i^k), then 
E[Ii] — 0. Hence the ICI power is, 



N-l 



(6) 



To appreciate the severity of frequency offset on the Carrier 
to Interference Ratio (CTIR) in an OFDM system, we define 
CTIR = £;[|J?;p]/jE[|/(p]. Assuming the channel transfer func- 
tion, H is unity and with zero AWGN, (i.e. Wi — 0), then 



CTIR = 



E^U' l/[iVsin7r(fc + € - Z)/iV]2 

[Nsinine/N)]^ 



(7) 



A similar analysis on the effect of frequency offset on OFDM 
can be found in [6] and [7]. Figure 2 shows CTIR in decibels 
as a function of e. Simulated values for AT = 64 and 256 using 
QPSK data mapping are also shown. It shows that the CTIR drops 
significandy when relative frequency offset, e increases beyond 
0.1, An important observation is that the CTIR is independent of 
A^. 

III. Schmidl and Cox Algorithm (SCA) 

One of the more robust schemes to estimate both frame syn- 
chronisation and frequency offset is the SCA. It uses two training 
symbols with the first one having a repetition within half a sym- 
bol period. Frame synchronisation is achieved by searching for a 
training symbol with two identical halves. If L = N/2, the sum 



of L consecutive correlations between pairs of samples spaced L 
apart is found as, 



L-l 



(8) 



n=0 



The output of the correlator, as given by (8) will reach a peak 
at the end of the CP of the first training symbol. If v > — I 
and p = V - (iV/i - 1), the peak is actually maintained for p 
samples just before the end of the CP of the first training symbol. 
Thus the correlator output will take the form of a plateau. This is 
because the last few samples of the CP, equal to the length of p, 
are not corrupted by the CIR. Timing synchronisation is achieved 
by locating the end of this plateau, denoted by dopt- Since the 
CP ensures that there is no Inter Block Interference (IBI) between 
OFDM symbols, the effect of the CIR in the first training symbol 
is cancelled if the conjugate of a sample from the first half is mul- 
tiplied by the corresponding sample from the second half. Hence 
the phase difference between the two halves of the first training 
symbol is caused by the frequency offset, 0 = ttc. It can be esti- 
mated Bis4> = l/p J2^=o ^P{(^opt - n) . If e < 1 there is no phase 
ambiguity in ^ and the frequency offset can be estimated as 



(9) 



In order to resolve potential ambiguity use can be made of the 
second training symbol, as detailed in [5]. The effect of AWGN 
on the partial estimation given by equation (9) is insignificant only 
if L, and hence the number of FFT points N, is large. The SCA 
performs well for OPTDM systems with in excess of 1000 [5]. 
For BFWA systems, data is transmitted in short bursts, particu- 
larly in the uplink. In this situation, it would be wasteful to use 
an OFDM system with large N. Besides large values of N will 
give rise to the additional problems of high peak-to-average power 
ratio and will also introduce latency which reduces protocol effi- 
ciency. Typically BFWA systems utilise lower values of A'' such 
as 64 to 256 as presented in [1] and [2]. When the SCA is used for 
systems with lower values of AT, the estimate e is not sufficiently 
accurate. This error results in a residual frequency offset that ro- 
tates the received constellation at a reduced rate, but one that is 
still significant enough to cause bit errors in coherently demodu- 
lated OFDM systems. In this paper we present the Residual Fre- 
quency Offset Correction Algorithm (RFOCA) that compliments 
the SCA by following it with a tracking function that continuously 
compensates for the residual frequency offset error. 

Other approaches that are used for frequency offset correction 
include for example the use of null symbols [8]. The use of a null 
symbol is not practical in a burst mode system as the transmission 
will always be a null while no data is transmitted. Other sugges- 
tions utilise the correlation between the cyclic prefix and the last 
few samples of the OFDM symbol [9] and the self-cancellation 
schemes presented in [10] and [7], where the same data is sent in 
more than one subchannel. Both schemes are not very practical 
when AT is a low value with the latter two approaches having the 
additional problem of reducing the number of useful subchannels 
by at least half. Note that the RFOCA does not utilise any pi- 
lots nor null symbols yet will be shown to effectively eliminate 
the residual offset. The proposal in [1 1] also attempts to track the 
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Fig. 3. Residual Frequency Offset Correction Algorithm (RFOCA) 

residual frequency offset. Here, the residual offset correction fac- 
tor is updated at the end of each OFDM symbol. However our tests 
indicate that the estimate is not very accurate when only one sym- 
bol is used, hence we update the residual offset correction factor 
only after analysing a block of A^^,, (A^^; > 1) OFDM symbols. 

IV. Residual Frequency Offset Correction 
Algorithm (RFOCA) 

In the proposed approach, an initial frequency offset acquisition 
is made using the SCA. We propose to estimate the residual fre- 
quency offset, e by tracking the rate of phase change at the FEQ 
outputs, Yrn.i ^s shown in Figure 3. However, the estimate of e 
could be seriously affected by subchannels with a low SNR result- 
ing from spectral nulls in the channel response Hi. An estimate 
of if; is made by taking the ratio of the transmitted and decoded 
output of the second training symbol of the SCA. Hence we only 
select those subchannels with \Hi \ above a certain threshold. We 
call this subset of subchannels d C [0, ..AT — 1]. The criteria ap- 
plied is to select subchannels with \Hi\ in excess of a standard 
deviation above the mean. For symbol m, the outputs of these 
subchannels y^^d are sent through a slicer to obtain Ym,d^ where 

ded. 

For symbol m the phase errors between Ym,d and Yrn,d^ namely 
Om,d are found as, 



0. 



m,d 



(10) 



The phase errors are stored in a buffer, 6 for a block of A^u, OFDM 
symbols. Note that only the selected subchannels will be stored. 
Hence if the number of subchannels selected in d is AT^ then the 
length required for the buffer is only N^N^. The Maximum Like- 
lihood (ML) estimate of the residual frequency offset e is calcu- 
lated as the gradient of the values stored in the buffer. Note, in 
this paper we have assumed QPSK data mapping though the tech- 
nique could be extended to other constellations. Initially, e can be 
quite high and if the values of the phase errors, O^^d > 7r/4, then 
it results in decoding errors being produced by the slicer which 
subsequently results in phase wrapping of 6m,d at ±7r/4. Hence 
it is imperative that we first unwrap Om^d before the calculation of 
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Fig.4. Phase Unwrapping Algorithm 

€. Although accurate phase unwrapping algorithnis are available, 
they are very complex. Since we are only interested in the phase 
gradient and not the exact phase values, a simpler scheme was 
selected for the unwrapping of the phase [12]. For purposes of 
clarity, we denote the ith sample of the wrapped phase and the un- 
wrapped phase as Oi and di, respectively. The phase unwrapping 
algorithm can be expressed as 



(11) 



where SAW(,) is a sawtooth function that limits the output to 
±7r/4 and a is a parameter that controls the variance of the un- 
wrapped phase. The algorithm is illustrated in Figure 4. 

The error due to frequency offset is compensated by multiply- 
ing the pre-FFT received symbols, r^.n with e^"^"*/^, where 
in general im is the frequency offset correction factor for sym- 
bol m, which is initialised to esc A at the acquisition stage. At 
symbol number qN^, iqN^ is updated according to ig^^ = 
iqN^-i -h iqN^, where q is an positive integer. For OFDM sym- 
bols [{q — l)Nyj,qNw — 1], i(q-i)N^ is used as the frequency 
offset correction factor. In other words im is updated once every 
iVt^ OFDM symbol blocks. 

The choice of is critical as a large value will cause the up- 
dating of ijn with too high a latency. However after a few updates 
the values of e are quite small, hence to get a better estimate in the 
presence of AWGN a larger value of will reduce the variance. 
The RFOCA can be summarised as follows. 

1. Do an initial frequency offset acquisition iscA using the SCA. 

2. Select the subset of subchannels with magnitudes exceeding 
the defined threshold in the channel transfer function, d, at the 
start of the burst. 

3. Obtain Ym,d from the FEQ output and then Y^.d by use of a 
slicer for each symbol m, where d £ d. Calculate and store Om,d 
for a block of iV^^, OFDM symbols. 

4. Find the unwrapped phase 0 from the wrapped phase 0, Find 
the gradient of 0, namely e. This process is done once every 
symbols. 

5. Once every A^^; OFDM symbol block, calculate the new fre- 
quency offset correction factor iq^^ . 

V. Simulation Parameters and Results 

OFDM systems with N = 64, 128, 256 have been simulated at 
a sampling rate of 20 MHz with a guard interval equal to 20 sam- 
ples, thus the subcarrier spacings are approximately 312 kHz, 156 
kHz and 78 kHz respectively. QPSK mapping for all subchannels 
has been employed and all the subchannels are used. A burst of 
320000 data bits is assumed to be transmitted. Each data point 
in the simulation results is obtained by averaging over 500 such 
bursts. In order to test the performance of the RFOCA with fre- 
quency offset alone, we have assumed perfect symbol and frame 
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synchronisation for the simulated coherent OFDM systems. A 
more robust symbol synchronisation algorithm that virtually guar- 
antees perfect frame synchronisation will be die subject of a fu- 
ture publication. After much testing, the values selected for N^v 
were 250, 250, 125 for N = 64, 128, 256, respectively. Figure 5 
shows the comparison of error variances at the end of the acquisi- 
tion stage using the SCA and the tracking stage using the RFOCA 
for A'' = 64, 128, 256 as a function of the Signal to Noise Ratio 
(SNR), with e = 0.5, It shows that the RFOCA error variance 
is many orders of magnitude lower than at the end of the acqui- 
sition stage. There appears to be a threshold effect in RFOCA 
at around an SNR of 16 dB, below which the error variance in- 
creases rapidly. It is mainly due to the errors present in the se- 
lected subchannels of d, which subsequently create estimation er- 
rors in the residual offset, Cm- Figure 6 shows the performance 
of the RFOCA in terms of Bit Error Rate (BER) vs SNR again in 
AWGN. The BER is virtually zero at values of SNR above 16 dB 
but suddenly increases below 16 dB. This result is a direct con- 
sequence of the threshold effect observed with the previous error 
variance results. It is observed that at an SNR of 16 dB only a few 
out of the total of 500 received bursts give rise to errors. How- 
ever, these errors are significant giving rise to a high overall BER. 
Such occurances became more common when the SNR is reduced 
below 15 dB, thus increasing the BER very rapidly. 

Appropriate models for BFWA channels are in the process of 
being defined. The Stanford University Interim (SUI) channels 
comprise 6 models for 3 different terrain conditions [13], All 
of them are simulated using 3 taps, each having either Rayleigh 
or Ricean amplitude distributions. The channel is assumed to be 
wide-sense stationary uncorrected scattering (WSSUS) and each 
tap of the CIR is modeled as hi = pie^^* , where the amplitude Pi 
and the phase (pi are selected independently [14], The burst takes 
less than 10 ms to transmit at the selected sampling frequency, 
consequently the channel is assumed constant for the duration of 
each burst. We have selected the SUI-2 channel model, pertaining 
to terrains with low tree densities and with antennas having a di- 
rectivity of 30 degrees at the SU and 120 degrees at the BS. The 
channel is characterised by a RMS delay spread of 0.2 //s. 

Figure 7 allows a comparison of the error variances of the SCA 
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and the RFOCA for an OFDM system with e = 0.5 in the SUI- 
2 channel. A new channel in accordance with the SUI-2 profile 
is randomly generated for the transmission of each burst. The 
RFOCA error variance is seen to be many orders of magnitude 
lower than at the end of the acquisition stage despite the presence 
of the SUI-2 channel. 

Figure 8 shows the performance of the RFOCA in terms of BER 
vs SNR for the SUI-2 channel and e = 0.5. The threshold effect 
now occurs at an SNR in the region of 20 dB owing to the effect 
of the SUI-2 channel. 

VI. Conclusion 

We have presented a Residual frequency Offset Correction Al- 
gorithm (RFOCA), that compliments the frequency acquisition 
process performed by the Schmidl and Cox Algorithm (SCA). The 
RFOCA continuously tracks and compensates for the residual fre- 
quency offset that is present after the acquisition stage. We have 
applied the RFOCA for AWGN and BFWA channels, namely the 
SUI-2 channel profile, and have shown a significant reduction in 
the error variance brought about by the use of the RFOCA fol- 
lowing the initial acquisition stage. Although the algorithm may 
appear to suffer from a threshold effect it still gives a significant 
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performance advantage at realistic signal to noise ratios. In the 
future we hope to address the threshold effect by improving the 
phase unwrapping and investigating other criteria for the selection 
of the subchannels. We also hope to improve the performance of 
the SCA in a frame synchronisation role. 
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